Quantum squeezing of optical dissipative structures 
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Abstract 

We show that any optical dissipative structure supported by degenerate optical parametric oscillators contains a special 
transverse mode that is free from quantum fluctuations when measured in a balanced homodyne detection experiment. The 
phenomenon is not critical as it is independent of the system parameters and, in particular, of the existence of bifurcations. This 
result is a consequence of the spatial symmetry breaking introduced by the dissipative structure. Effects that could degrade 
the squeezing level are considered. 



PACS numbers: 42.50.Lc, 42.65.Sf 

Introduction. Vacuum quantum fluctuations consti- 
tute the ultimate noise source affecting any coherent ra- 
diator, like a laser. These fluctuations define the so-called 
standard quantum limit as they set the maximum pre- 
cision attainable with classical optical techniques, even 
rendering the latter useless in some applications such as 
precision metrology Q and quantum information proto- 
cols |2|, y] ■ It is possible however to break this limit with 
the help of quantum states of light: Squeezed states , 
displaying fluctuations below the standard quantum limit 
in one of the field quadratures, play a prominent role in 
this regard, and are by now routinely generated, e.g., 
by single- mode optical parametric oscillators/ amplifiers 

Bala. 

In the last fifteen years a new branch of quantum optics 
has emerged under the general name of " quantum imag- 
ing" that studies the spatial aspects of the field quan- 
tum fluctuations and generalizes the study of squeezing 
to multimode beams [5j- We consider in this context 
the squeezing properties of optical dissipative structures 
(DS) -transverse patterns-, which are self-sustained, sta- 
ble spatial structures that form across the plane perpen- 
dicular to the axis of multi transverse-mode nonlinear 
resonators 0. We consider two outstanding classes of 
optical DS in the planar degenerate optical parametric 
oscillator (DOPO): periodic patterns and localised struc- 
tures. 

Localised structures of nonlinear optical resonators are 
also called cavity solitons (CS) for their resemblance with 
optical fibre solitons and with optical spatial solitons. 
However, in spite of their appealing resemblances, the lat- 
ter are Hamiltonian objects resulting from perfect com- 
pensation between dispersion/diffraction and nonlinear- 
ity, while CS are also ruled by dissipation and feedback 
and are not true solitons in the mathematical sense. The 
squeezing properties of optical fibre solitons are well un- 
derstood since long time ago 0, E| , and those of optical 
spatial solitons have been considered more recently [lfT |. 

We show below that any stationary DOPO DS is 



phase-squeezed in a special transverse mode in the linear 
approximation. Importantly, the degree of squeezing is 
independent of the system parameter values and of the 
existence of bifurcation points. This result is a direct 
consequence of the free diffusion of the DS across the 
transverse plane, ruled by quantum noise. Although dif- 
ferent in nature, the reported phenomenon resembles the 
quantum noise suppresion on the difference of intensities 
(amplitude-squeezing) of a two-mode optical parametric 
oscillator above threshold , which is associated to the 
existence of a continuous diffusion of the phase difference 
between the two modes. In our case the fact that squeez- 
ing occurs in a transverse mode with a special spatial 
shape could be useful for some applications Q. 



Model. We consider the model of for a DOPO 
with plane cavity mirrors. A plane wave coherent field 
of frequency 2u> s and amplitude £;„ pumps the resonator 
containing a crystal, which converts pump photons 
into signal photons (of frequency lo s ) and vice versa. 
Only two longitudinal cavity modes, of frequencies loq 
(pump mode) and u>\ (signal mode), the closest to 2lu s 
and u> s , respectively, are assumed to be relevant. These 
modes are damped at rates 7„ (n = 0, 1) and losses 
are assumed to occur at a single cavity mirror. The 
intracavity field envelope operators for pump and sig- 
nal modes are denoted by Aq (r,t) and A\ (r, t), respec- 
tively, where r = {x, y) denotes the transverse coordi- 
nates, which obey standard equal-time commutation re- 
lations [A n (r, t) , A n fr', t)] = 5 (r — r'). Using the pos- 
itive P- representation IT3L fl4| . which sets a corre- 
spondence between the quantum operators A n and A n 
and the independent, stochastic c-number fields A n and 
A„ , respectively, we derive the following model equations 
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have the form: 



d t Ao (r, t) = -70 (1 + iA ) A + £ in + iy'?V 2 A (1) 
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d t Ai (r, t) = -7i (1 + iAi) A + 271^ V 2 A 



(2) 



+ gA^[ + r? (r, t) , 



where Ao = [ojq — 2w s ) /70 and Ai = (wi — u) s ) /71 
are cavity detunings, /1 = c/sflux^jl is a characteristic 
(diffraction) length, c is the speed of light in the crystal, 
V 2 = d 2 /dx 2 + d 2 /dy 2 accounts for diffraction, and g is 
the (real) coupling coefficient proportional to the relevant 
second-order nonlinear susceptibility of the crystal. As 
everywhere along the rest of this Letter the equations for 
A n are to be complemented by those for the "hermitian 
conjugate" fields A+, which arc obtained from those for 
A n by complex-conjugating the parameters and by doing 
the replacements A n < — ► A^ and 77 < — > rj + . Finally r\ 
and rj + are independent, real white Gaussian noises of 
zero mean and correlations 



(7/(r,i)77(iV)) 



(r, + (r,t)r]+(v',t')) (3) 
5 (r - r') 5(t- t') . 



Thus A n and A„ are not complex conjugate, al- 
though the stochastic average of any function of them, 
(/ (Ao,Aq,Ai,Ai)), yields the normally ordered quan- 
tum expectation value ^: f (^Ao, Aq, A\, A\\ :Y For ex- 
ample, the signal photon density is calculated as 

M (r, t) = (A\ (r, t) A, (r, t)) = (A+ (r, t) A 1 (r, <)) . 

(4) 

In the limit of large pump detuning (70 |Ao| ^> 
71 I Ai| , 70, 71) pump diffraction plays a negligible role 
and pump fields can be adiabatically eliminated [l6| 
as Ao = 71 (fi + i-jpAf) /g, where n = \A \/g, 
fi = g I fi„ I / (71 I Ao|) > is the dimcnsionlcss pump pa- 
rameter, a = sign Ao, and we wrote £ m = ia \S- ln \ without 
loss of generality. Then eq. (J2J) becomes 



9Mi=7i - (1 + *Ai) Ai + fxAf 



(5) 
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Classical dissipative structures. When noises are 
ignored and Af is identified with A* (classical limit) eq. 
(0) coincides with that for a classical DOPO with large 
pump detuning [J^, which supports differen t ty pes of 
steady and stable DS, both periodic patterns |l7| or CS 
[ill fig depending on the parameter region. All these DS 



Ai(r)=Ai(r-n), A (r) = ree lCT V (r) 
e 2la9 = M " X (l + J( tv / M J— l) , 
(al\ V 2 - f + T 2 ) T = 0, [3 2 = a At + ^fW~- 



(6) 
(7) 

(8) 



where ri = {x\,y\) is arbitrary due to the translation 
invariance and T is real. Note that the only parameters 
defining the DS are {er, /1, Ai} as re and l\ merely act as 
scale factors. 

Dynamics of quantum fluctuations. Fluctuations 
around any classical DS, eq. @, are studied by setting 

A 1 (r,t)=A 1 (v-r 1 ) + a 1 (r-r 1 ,t), (9) 
At(r,t) = A(r-ri) + <(r-ri,t), (10) 

where, given the translation invariance, the position of 
the classical DS, ri (t), is let to vary in time as it is an 
undamped variable that is excited by noise. Expressing 
the fluctuations as a = (<zi,a^) (T denotes transpo- 
sition) and linearizing the resulting eq. stochastic 
equations for the quantum fluctuations are obtained: 



— re (G x dxi/dt + G y dyi/dt)+d t a = 71 Ca- 
where 



Tlh, (If) 



h = I V«o V, 



Gj = 



'd 3 T, j = x,y, (12) 
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-1 - iLi 



a -- 

-1-iLi 

d 
-1 + iLi 



■ 72 



(13) 
(14) 

(15) 



with L\ = (l 2 V 2 — Ai). The spectra of the linear op- 
erators C and £t, which we introduce as £vj = AjVj, 
£t Wj — A*Wi, are clearly relevant 0. In order to deal 
with stable DS, it is assumed that ReAi < for any i 
. With the usual definition of scalar product (b| c) = 
J d 2 r (r) • c(r), the relation (wj|£c) = A,(wj|c) 
holds. We assume that all eigenvectors are suitably or- 
thonormalised as (wj|vj) = Sij. In general the spec- 
tra must be computed numerically; nevertheless two gen- 
eral properties of the discrete spectra can be stated: (i) 
G x ty-\ are Goldstone modes as CG x r y ) = (we write 
Vi x ily) = G x ( y ) and denote by w lx ( l9 j the associated 
adjoint eigenvectors: /^Wwi^) = 0); and (ii) 



-2w 2a ;(2j/)) w 2x(2y) — (w 2x ( 



(2y), w 2x{2y) 



W2x(2y) 



(y)- 



(16) 



Properties (i) and (ii) occur independently of the set 
of parameters {cr, /x, Ai}. Property (i) is a mere con- 
sequence of the translational invariance of the problem. 
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Property (ii) is the key for our analysis. Note that 
these eigenvectors exist as the classical DS breaks the 
spatial symmetry: for a spatially homogeneous solution 
(T = est) G x ( y} = 0. 

Before solving eq. fTTt we consider its projections onto 
the eigenvectors Wwi y ) and W2x(2 V ) introduced above: 



dx\/dt 



7i K 1 £lx> 



dc 2x /dt = -27x020, + V7i^2 



(17) 

where & (t) = (Wj| h) are noise sources, C2 X (t) — (w 2x \ a) 
is the projection of the fluctuations onto the eigen- 
mode W2, and corresponding expressions for dyi/dt and 
dc 2y /dt. Note that the equation for x\ is diffusive, as an- 
ticipated, because Goldstone modes are excited without 
cost [i^- We further notice that the equation for ci x is 
analogous to that derived in [2l| for the hexagonal mode 
stationary phase in a Kerr cavity, which is later inter- 
preted as the hexagonal pattern transverse momentum 
in H2. 

Squeezing via optical homodyning. We consider 
the squeezing properties of a DS as measured in a bal- 
anced homodyne detection experiment (2^: The outgo- 
ing quantum field, Ai jOUt (r, t), is combined in a beam 
splitter with a local oscillator field (LOF) that lies in an 
intense (multimode) coherent state of transverse complex 
envelope ql (r — i"l (t)), which is allowed to be dynami- 
cally shifted. In the detection of squeezing one measures 
the normally ordered part of the fluctuation spectrum of 
the intensity difference between the two output por ts of 
the beam splitter, S, which can be computed as [23j 

S H = 2 7 i / dTe- iuT {6£ H (t + t) S£ r (<)), (18) 



S£k (*) = 



J d 2 r |a L | 



:<a L (r + p(t))|a(r,t)), (19) 



a h = {a h ,a L ) , p = n - r L . 



(20) 



When A ljOUt (r,i) is in a (multimode) coherent state 
S (lu) = 0, the standard quantum limit. On the other 
hand S (u s ) = — 1 signals complete absence of quantum 
fluctuations at to — u> s . 

Let us assume momentarily that we can set p = 
in eq. (|19|) . which means that we can shift the LOF 
according to the DS movement. Let us choose a LOF 
with ql = i>G x (y) (i-e., ol~l =2x, eq- (HD) so that 
<5£h (t) = c 2x (i), see after eq. (|T7|) . Standard techniques 
applied to eq. (|17fl allow to compute the stochastic 



13] 



correlation (<5f H + r ) (*)) = — ^ ( 
using eq. (fTHf> we get 



-2 7 l U 



47? 



Finally 



(21) 



which is the main result of this Letter: As S (oj = 0) = 
— 1, DOPO DS display perfect squeezing at u> = when 
probed with the appropriate LOF (q;l = iG x ^). As eq. 



(|21|l is independent of the kind of DS and of the system 
parameters, the result is universal and independent of 
the existence of bifurcations. This LOF is, in principle, 
easily realisable as it is the 7r/2 phase-shifted gradient of 
the corresponding DS envelope, eq. ©, which can be 
easily synthesised by, e.g., Fourier filtering. 

It is interesting to notice that in |2l|. a perfectly 
squeezed spatial mode was identified in the hexagonal 
pattern arising in a Kerr cavity and, as in our case, the 
result is independent of the parameter values. Although 
derived by different means from the ones used here, this 
result is very likely connected to the one we have just 
derived. 

We note that the linearised approach is valid, in prin- 
ciple, when all eigenvalues are strictly negative, as then 
all fluctuations remain small. In our case however a null 
eigenvalue exists always -that associated with the Gold- 
stone mode. Nevertheless that eigenvalue is just the re- 
sponsible for the continuous diffusion of the position of 
the DS (similarly to the continuous diffusion of the phase 
difference in |ll|), and does not entail an energetic diver- 
gence. Hence one can be confident that the linearised 
theory developed here represents quite an accurate de- 
scription and that a nonlinear treatment jl5j would not 
lead to dramatically different results. 

Next we consider two effects that could degrade the 
measured squeezing level. Although everything to be said 
applies to any DOPO DS, we focus on the bright CS, 
which exists for a = +1 1G, 18], for the sake of clarity. 

Influence of the CS movement. Equation l|21() 
is valid if we use a movable LOF which exactly follows 
the CS movement. This could be done by tracking the 
movement of Ao j0 ut, which is correlated with j4i )OU t, with- 
out disturbing the subharmonic CS. The output of this 
continuous measurement would be then fed into a posi- 
tioning system controlling i"l (t), giving rise, in general, 
to a time delay so that tl (t) — ri (t — td) yielding 
pit) — V\ (t) — vx (t — td). The point is how much the 
CS position diffuses in time as compared with its width 
Ax. Standard techniques 0] applied to eq. lfT7|) al- 
low to obtain (p 2 (t)) = Dt^, where the diffusion con- 
stant D - 7i k" 2 Q. Thus (p 2 (i)) /Ax 2 ~ Jit d /N u 
being Ni ~ (kIi) 2 the number of intracavity signal pho- 
tons in one CS and Ax ~ l\ its width 25]. We see 
that Nx acts as an inertial mass. Using realistic values 
for the system parameters [2(| one has N\ ~ 10 12 , and 
v/(p 2 {t))/Ax < 5 • 1CT 4 for a delay time i d < 1ms. 
Thus the relative error existing between the location of 
the CS center and that of the LOF is very small as 
compared with the CS width. In order to assess the 
negligible influence of this effect we expand eq. I|19(l 
up to second order in p. The squeezing spectrum is 
then given by eq. I|21|) plus a correction proportional 
to (p 2 (t)) I {Axf < 3 ■ 1CT 7 , again for t d < 1ms, which 
is absolutely negligible. 
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FIG. 1: Squeezing level (at the labeled frequencies) displayed 
by the ID CS when nonideal LOFs are used. In (a) a Gauss- 
Hermite LOF (GHi) of width £ is used. In (b) LOFs displaced 
x from the CS center are considered (w2 denotes a special 
LOF, see text). Ax = hf/3 denotes the CS width. Parameters 
are a = +1, Ai = 1, fx = 1.2. 



Let us now consider how the squeezing properties of 
the CS are modified if the LOF is kept fixed, which cor- 
responds to a simpler scheme. Due to the unbounded 
movement of the soliton, one must perform the hetero- 
dyning experiment in a short time (call it tu) in order to 
obtain significant squeezing. If we assume that at t = 
the LOF and the CS centers are made to coincide one has 

P(t) 
lfj- 6 - 



= ri (t) - n (0) and y/(^(tj)/Ax ~ ^pi ~ 
Ayii. Then, if we take tn < 1ms (for the used pa- 
rameters) we can again take p = in eq. (|19|1 . which 



yields (S£ R (t + r) 5£u (*)) 



i e -2 7 i|r| as before. The 



squeezing spectrum is then given by eq. I|18(l with the 
limits of integration being replaced by T- 2 • The result 
reads as eq. (|21|l plus a correction proportional to e~ 7l * H , 
which is virtually zero. Then the obtention of almost op- 
timal levels of squeezing is not affected in practice by 
the existing CS movement. We note that this insensi- 
tivity contrasts with the issue of quantum images |27| . 

whose squeezing properties are washed out by their jitter- 

l li 

ing. This is a consequence of the strong inertia (oc N 1 ) 
that the CS movement displays against fluctuations, as 
compared with the below- or close-to-threshold emission 
analysed in [27|. 

Influence of the shape and positioning of the 

LOF. Up to here we have dealt with a special LOF. 
However the use of a LOF with the exact form or pre- 
cisest positioning is not critical as we show next. In 
what follows we ignore the negligible influence of the 
CS movement and consider a static LOF. The study re- 
quires fully solving eq. I|ll(l . what was done by using 
the (biorthogonal) basis {vi,Wj} formed by the eigen- 
vectors of the linear operators L and & in eq. Q14JI 28] . 
This technique is highly convenient as it allows to circum- 
vent the numerical simulation of the Langevin eq. (|llfl , 
which can be problematic and, in any case, extremely 
time consuming. We write the field fluctuations and the 
LOF vector as a(r, t) — Yl°i (*) v » ( r )> °i = ( w il a )> 
where the expansion excludes the Goldstone modes as 
ri is to denote the position of the CS. By substituting 



this expansion into eq. (|llf> and projecting, one obtains 
dci/dt = -fiXiCi + ^/tT (wj| h). The study is further fa- 
cilitated by expressing a general LOF as c*l = J2 a i w ii 
oti = ( v i| <*l)- One obtains 



S(u>) 

Sp,q H 



27i 



/ d 2 r |a L | p,<? 



dre 



^2 a * P a i s p-i H ' ( 22 ) 

C p (t + T) C q (t)) , 



which are easily evaluated by using standard methods 
[l3j |. In the end all we need is to compute the spec- 
tra of L and C\ which was done by using a Fourier 
method [2^. We limit here our study to the ID CS 

for the sake of computa- 
tional economy. The influence of the LOF shape was 
studied for a Gauss-Hermite mode of appropriate phase, 
t*L {x) — GH\ (x) = ie t6 x e~2W0 , which is simi- 
lar to W2 — iG x — ie i d x T. Figure 1(a) shows that 
quite high levels of squeezing can be reached even with 
this non ideal LOF. Finally the influence of a misposi- 
tioning was studied both for ql (x) = W2 {x — x) and 
(*l {x) = GH\ (x — x), fig. 1(b): Mispositionings as high 
as 15% of the CS width still yield quite good levels of 
squeezing as well. 

In summary we have shown that optical dissipative 
structures sustained by a DOPO always contain a trans- 
verse mode that is completely free from zero-frequency 
quantum fluctuations. Unlike single-mode cavity squeez- 
ing, which is perfect only at bifurcation points, our result 
does not depend on the parameter setting. 
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